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ABSTRACT 

We present a theoretical model fit to the ffi>T/STIS optical transit transmission spectrum of 
HD209458b. In our fit, we use the sodium absorption line profile along with the Rayleigh scat- 
tering by H2 to help determine the average temperature-pressure profile at the planetary terminator, 
and infer the abundances of atomic and molecular species. The observed sodium line profile spans an 
altitude range of ^3,500 km, corresponding to pressures between ^0.001 and 50 mbar in our atmo- 
spheric models. We find that the sodium line profile requires either condensation into sodium sulfide 
or ionization, necessary to deplete atomic sodium only at high altitudes below pressures of ~3 mbar. 
The depletion of sodium is supported by an observed sudden abundance change, from 2 times solar 
abundance in the lower atmosphere to 0.2 solar or lower in the upper atmosphere. Our findings also 
indicate the presence of a hot atmosphere near stratospheric altitudes corresponding to pressures of 
33 mbar, consistent with that of the observed dayside temperature inversion. In addition, we find a 
separate higher altitude temperature rise is necessary within both the condensation and ionization 
models, corresponding to pressures below ^0.01 mbar. This hot higher altitude temperature indi- 
cates that absorption by atomic sodium can potentially probe the bottom of the thermosphere, and 
is possibly sensitive to the temperature rise linked to atmospheric escape. 

Subject headings: planetary systems - radiative transfer - stars: individual(HD209458) 



1. INTRODUCTION 

Early theoretical tra nsit models of hot-Jupiter exo- 
planetarv atmospheres (ISeager fc Sasselovl [20001 : iBrownl 
12001b iHubbard et al.ll2001l ) contained strong sodium ab- 
sorption signatures, predicted to be as high as 0.2% for 
clear cloudless atmospheres. The original measured Na 
absorption signat ure in HD209458b f r om m edium reso- 
lution HST data ICharbonneau et all (|2002l ) was, how- 
ever 1 _j^><smBJler than initial predictions. Addition- 
ally, (Barman] (|2007f ) has also identified a potential broad 
band Na signature from photometry of low resolution 
STIS data in the sam e planet, taken several years later 
(|Knutson et al.l l2007f h These observations have lead 
to multiple theories to explain the small Na signa- 
ture. The most widely cited explanations of reduced 
Na absorption include: high clouds or hazes, ionization, 
low global Na abundance, ra inout, as well as non-local 
thermal equilibrium effects (Charbonncau ct al 120021 : 
iBarman et aTll2002t iBarmad 120071 : IFortnev et al.ll2003f ). 
High opaque clouds could potentially mask absorption 
features by effectively cutting transit information for al- 
titudes lower than the c loud deck (|Seager &: Sasselovl 
120001 IFortnev et ai1l2003h . The small atmospheric sig- 
nature has made detection from ground-based telescopes 
difficult, with no Na detection cur rently reported for 
HD209458b after m ultiple attempts (|Narita et alj 120051 : 
lArribas et~aill2006f) . 

In the grazing line of sight of transit geometry, the mea- 
sured effective altitude depends upon the wavelength de- 
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pendant opacities of the atmospheric constituents giving 
a transit transmission sp ectrum its altitude sens itivity 
(see relevant equations in iLecavelier et all 1200 8 a). The 
effective altitude is also dependant upon the abundance 
of absorbing atmospheric particles, as well as the atmo- 
spheric pressure. As the abundances are not in general 
a priori known, this leads to atmospheric models which 
have degenerate solutions between modeled abundance 
and pressure, where an atmospheric transit model with 
increased abundance is identical to models with lower 
abundances and higher pressures. 

Here we interpret the optical HST STIS transmission 
spectra of HD209458b of lSing et all (|2OO80 . The spectra 
consists of STIS data obtained during planetary tran- 
sit at low and medium spectral resolutions. The two 
datasets are combined to extend the measurements over 
the entire optical regime, providing a way to simulta- 
neously measure, in the case of atomic Na, both the 
narrow line core and wide line wing absorption as well 
as other absorbing or scattering species. Since the Na 
line core and wings are optically thick at different at- 
mospheric heights, a detailed analysis of the Na line 
profile in conjunction with other observed species can 
potentially break the model degeneracies and reveal at- 
mospheric structure such as the altitude-temperature- 
pressure (z-T-P) relation, abundance, and the possible 
presence and height of planetary clouds. In this paper, 
we focus mainly on the resulting temperature pressure 
profile, as well as the Na line itself. We briefly describe 
the observations and resulting spectra in §2, present our 
transit model fits in §3, report on an observed tempera- 
ture inversion in §3.3, and dis cuss our results in§4_& 
5. In addition to this work, ILecavelier et alj (|2008bl ) 
provides a detailed de s cripti on of H2 Rayleigh scatter- 
ing while iDesert et all (|2008l ) goes into modeling details 
of TiO, VO, and other molecular species. 
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2. OBSERVATIONS 

The data reduction, procedures, and various tests re- 
sulting in the transmission spectra use d in th is analy- 
sis are covered in detail by iSing et a"T1 (|2008| ), though 
we briefly describe them here for completeness. The 
ffST/STIS optical transit data for HD209458b consist 
of data at low resolution covering ~3, 000-10,0001 with 
the G430L and G750L gratings and medium resolution 
data from the G750M grating covering 5,813-6,3821. 

With the medium resolution data. ISing et al.1 ((2008) 
measured the Na absorptio n using the double differen - 
tial photometric method of iCharbonneau et al.l (|2002f) . 
as well as a double differential spectroscopic evaluation. 
For th e differential photometric measurement, Sing et al. 
( 2008) selected 23 wavel ength bands from 4 t o 100 A cen- 
tered on the Na D lines, ICharbonneau et all (|2002D used 
three different wavelength bands, and compared the ab- 
sorption to that of regions red-ward and blue- ward of the 
Na D line region. The resulting photometric Na absorp- 
tion profile (see Sing et al. 2008, Fi g 3) matches the 
original results of ICharbonneau et al.l (|2002f ) and, in ad- 
dition, reveals the strong absorption from the Na D line 
centers. 

To analyz e the low and med ium resolution HST spec- 
tra together. ISing et ail (|2008l ) used a limb-darkened cor- 
rected transit spectral ratio to probe the full wavelength 
dependence of atmospheric absorption features (see Fig. 
1 & Fig. 2). From this analysis, the Na D wavelength 
region shows an absorption plateau between ~5600-6700 
A, which is ~0.045% above the minimum absorption level 
observed (see Fig. 1). Note that the wavelength range of 
the medium re solution data, originally u sed to detect at- 
mospheric Na ()Charbonneau et a.l.ll2Q02h , is entirely con- 
tained within this broad plateau absorption region. The 
total absorption seen in the Na D line cores, therefore, 
is the result from the medium resolution data, which is 
the absorption above the plateau region measured to be 
^0.065%, plus the absorption level of the plateau region 
itself. This total Na core absorption level is ~0.11%, 
within a 4.4 A band, and is much closer to the o riginal 
predictions (|Seager fc Sasselov|[2000HBrownll2001h . indi- 
cating that a large extent of the atmosphere is being 
probed during transit. Additionally, the Na D doublet is 
resolved in our resulting spectra which shows two strong 
narrow absorption lines (see Fig. 2). 

2.1. Observed Na line profile 

The observed Na D transit line profile (Fig. 1), 
with it's narrow absorption core and broader plateau- like 
shape, is in stark contrast to cloudless theoretical mod- 
els of HD209458b (fSeager fc Sasselovl l2000l: lBrownll2001l: 
iFortnev et al.|[2003h lBarmanll2007h which have smoothly 
decaying Lorentzian "Eiffel Tower" shaped line wings. 
This more complicated observed line profile suggests a 
very low quantity of Na is present in the middle at- 
mosphere and much larger quantity in the lower atmo- 
sphere, pointing to either high altitude Na condensation 
into Na2S or alternatively, depletion via ionization. 

Previous studies have predicted models of HD209458b 
which included Na 2 S clouds (jlro. Bezard. fc Guillotl 
2005), as high altitude temperatures on the night side 
reach levels cool enough for condensation of atomic Na 
(|Lodderdll99lH Sh owman fc Guillotl 12001) . Sodium sul- 



fide would effectively decrease the amount of atomic 
sodium at altitudes above the condensation tempera- 
ture, creating a stratified-like atmospheric profile de- 
pleted of atomic Na only at those high altitudes. The 
transit signature with such a layered Na profile produces 
a broad plateau shaped line wing along with narrow ab- 
sorption cores rising above t he plateau to higher levels 
(jlro. Bezard. fc Guillotl [20051) . Na 2 S would have to be 
transparent, since they are not observed in absorption 
during transit at the limb. If Na 2 S or other clouds/hazes 
were opaque in absorption during transit at the conden- 
sation altitude, significant absor ption would ex tend over 
a very broad wavelength range ((Fortnevi 12005). masking 
the low resolution absorption features we observe. Such a 
flat, largely featureless transmission spectra has recently 
been observed in HD189733b, the result of hig h altitude 
haze (|Pont et al.ll2008t iLecavelier et al.ll2008ali . 

An alternative metho d to deplete neutr al atomic 
sodium is by ionization (|Fortnev et al.l 120031 ). In this 
scenario, UV photons between wavelengths of 1,107 and 
2,412 A (corresponding to 11.2 and 5.14 eV respectively) 
photoionizc sodium along the tra nsit line of sight unti l 
recombination matches ionization. IFortnev et al.l (|2003l ) 
found a sharp boundary at the limb between higher at- 
mospheric regions of complete ionization and lower re- 
gions with no ionization. Such a model could explain the 
observed Na line shape, with a narrow Na absorption core 
rising above broad line wings, if the atmosphere suddenly 
becomes ionized above a certain altitude. 

3. MODEL FITS TO THE TRANSIT TRANSMISSION 
SPECTRA OF HD209458B 

We modeled the atmospheric transmission spectra cal- 
culatin g the transit signatur e based on a theoretical 
model (|Ehrenreich et al1l2006l ) containing absorption by 
Na, H 2 Rayleigh scattering, TiO and VO. The model out- 
put absorption spectrum can then be directly compared 
to our limb-darkened corrected transit absorption spec- 
trum. We calculated the optical depth along a grazing 
line of sight through the planet's atmosphere, and calcu- 
lated the absorpt ion depth at a g i ven wa velength using 
equations given in lLecavelier et all ((2008a) . For the Na D 
lines, we used collis ion-broadened line-shapes cal culated 
using Equations of llro. Bezard. fc Guillotl (|2005[ K Con- 
cerning atomic sodium, ultimately the free parameters 
to fit in the transit calculation include the atmospheric 
temperature profile as a function of altitude, and the 
mixing ratios of Na below and above the condensation 
or ionization level (see Fig. 3). 

We found two families of solutions, with each giving 
similarly good fits to the data. In one model, Na is 
allowed to condense into sodium sulfide at the chemi- 
cal equilibrium temperature, while in the other model 
neutral Na is depleted at high altitudes, presumably by 
ionization, below an arbitrary pressure. 

3.1. Condensation model 
3.1.1. Temperature- Pressure profile 

For the calculation, we fit various simplified T-P pro- 
files in which the temperature is taken as a linear func- 
tion of the altitude and pressure is calculated using hy- 
drostatic equilibrium. Note that it is not specifically the 
T-P profile shape itself which is important in our fitting 
of the data, but rather the different temperature and 
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pressure regimes imposed by the observations. At our 
signal-to-noise levels, the general variation of the T-P 
profile in between the different atmospheric regimes we 
find (hot at high pressure, cool at middle pressures, and 
hot at low pressures) is only loosely constrained. A more 
complicated T-P profile would add undesired additional 
parameters into the fit, leading us to choose as simple a 
T-P profile as possible. 

We found it necessary to have a minimum of three 
variable points in the T-P profile. The temperature as a 
function of the altitude, T(z), is thus composed of two 
linear functions below and above the altitude z m , where 
z rn is the altitude of the minimum temperature, T m , We 
took the lower reference altitude, z s , for t he altitude cor- 
respo nding to the mean planetary radius (|Knutson et al.l 
l2007fl . and the altitude z t h is the fixed position, at high 
altitude, needed to define the second part of the T(z) lin- 
ear function. We used z t h corresponding to the altitude 
of the absorption depth of 1.55%, which corresponds to 
the mean absorption depth in the core of the sodium 
doublet. The temperature and pressure at the altitude 
z s , the altitude z m of minimum temperature T m , and 
Tth are the five free parameters defining the T-P pro- 
file. The sodium abundance in the lower and the upper 
atmosphere are the two other free parameters in the fit. 

In our condensation fit we found: (1) a high tempera- 
ture of 2,200±260 K at 33±5 mbar (~ 500 km) mainly 
constrained by Rayleigh scattering (|Lecavelier et al.l 
2008b) seen at wavelengths shorter than ^5,000 A, (2) a 
low temperature of 420±190 K at 0.63±0.27 mbar needed 
to explain the plateau at 0.045% (~1,500 km) above the 
minimum absorption level as a Na abundance change due 
to condensation, and finally (3) a high altitude (~3,500 
km) high temperature of 1,770±^ K at 0.0063±0.003 
mbar needed to fit the strong peak and narrow width of 
the cores of the Na I doublet. 

The pressure-temperature gradient we find when im- 
posing condensation is close to and consistent (within 
la) of an adiabatic gradient. A temperature gradient 
which is adiabatic or less is needed to be stable against 
convection, and solutions imposing adiabatic variations 
give acceptable fits to the data (an adiabatic fit is shown 
in Fig. 3). We ultimately, however, choose to quote val- 
ues found from the more general solution, as the linear 
temperature-altitude relation we assume includes adia- 
batic gradients and the general fit is located at a well 
behaved minimum y 2 . Although the best fit solution is 
slightly superadiabatic, the relatively low quality of the 
data still allows for a wide range of T-P gradients and the 
adiabatic-like profile found from our general fit indicates 
that an adiabatic behavior is likely at those altitudes. 

3.1.2. Atmospheric Na Condensation, TiO & VO 

Our model contained two different amounts of atomic 
Na having separate variable parameters for the mixing 
ratios of the middle and lower atmospheres. The pressure 
boundaries (or condensation pressure) where the quan- 
tity changes occur is taken to be at an altitude where the 
temperature profile crosses the Na2S condensation curve. 
In our fits, we found a Na mixing ratio of 3.5j^'g x 10 - ' 

in the middle atmosphere and 4.3^'g x 10~ 6 in the lower 
atmosphere, corresponding to 0.2 and 2 times solar abun- 
dance. 



IFortnevI l|2005h estimated that the optical depth of 
Na2S condensates in transit geometry could be optically 
thick for HD209458b (r ^0.73) for an atmosphere of so- 
lar composition and cloud base at 30 mbar. Given the 
fit T-P profile, we find a Na2S cloud base of ~3 mbar. 
Lower pressures or larger particle sizes would provide a 
natural explanation for transparent clouds, as the con- 
densate optical depth is proportional to th e pressure an d 
inversely proportional to the particle size (|Marleyll2000l ). 
Given our determined Na abundances, the mixing ratio 
of Na2S we find is (2±l ) xl0~ 6 , which matches the as- 
sumed value of Fortncy 1 20051) . Scaling the estimated 
optical depth of IFortnevI (|2005h . a cloud base at 3 mbar 
gives r ~0.073, indicating that transparent Na2S con- 
densates are plausible. 

The continuum blue-ward of the Na D lines, between 
~4,000 and 5,500 A, shows a slightly lower absorption 
than the red-ward region, 7,000-8,000 A, and a signifi- 
cant absorption 'bump' around ^6,250A which appears 
in both resolutions. Our fits indicate that the red-ward 
absorption levels (7,000-8,000 A) are likely due to TiO 
and VO at stratospheric heights, being largely confined 
to altitudes below our Na condensation level (also see 
iDesert et all 120081) . The TiO and VO features mask 
part of the long wavelength Na line wing, though not 
the short wavelength side, providing strong constraints 
on the abundance and altitude of TiO and VO. The hot 
low altitude temperatures we find from Rayleigh scatter- 
ing are consistent with the presence of TiO and VO, as 
both species begin to condense out of the atmosphere at 
temperatures below ~ 1,800 K. 

The low-altitude pressure assumed h ere is dependent 
on the detection of Rayleigh scattering (jLecavelier et al.l 
2008b) , though multiple interpre tations of the absorption 
rise at those wavelengths exist (fLecavelier et al.l [2008b; 
iBallester. Sing, fc Herbert! l2007t iBarmanl |2007| ). Ulti- 
mately, it may be difficult without Rayleigh scattering to 
assign an absolute pressure scale to these measurements 
which are inherently differential. The effective planetary 
radius seen during transit is sensitive to the quantity 
of the atmospheric pressure tim es the abundance of ab - 
sorbing species (see equations of lLecavelier et al]|2008al) , 
making an atmospheric fit at higher altitude possible by 
simply increasing the absorbers abundance. The advan- 
tage of detecting Rayleigh scattering by H2 when in- 
terpreting these models is that the abundance of H2 is 
known, constituting the bulk of the atmosphere, which 
thereby fixes the pressure, and defines the pressure scale 
for the rest of the transmission spectrum. Our adoption 
of Rayleigh scattering here is based on the goodness of 
fit to the data and the necessity to probe a large pressure 
range, given the total Na profile observed. 

Since the low-altitude pressure is an upper limit con- 
strained by Rayleigh scattering, a lower pressure would 
have an effect on the lower temperature profile, increase 
the lower Na mixing ratio, and decrease the condensate 
pressure determined in our transit fit. Whatever the true 
pressure P s or additional absorption sources, however, 
the basic atmospheric characteristics of Na determined 
here would remain the same. Namely, the middle at- 
mosphere has a greatly reduced amount of Na with a 
temperature below Na2S condensation, while the lower 
and higher atmospheres are warmer than the condensa- 
tion temperature with the lower atmosphere containing 
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a solar-like Na abundance and higher one a significantly 
reduced Na abundance. 

3.2. Photoionization model 

Given the large uncertainties modeling hot- Jupiter 
Na photoionization, we choose not to employ a com- 
plex ionization mod el, but rather follow the results of 
iFortnev et alj (|2003ft and assume in our model that neu- 
tral Na can be suddenly ionized above a given altitude 
leading to the observed Na deplet ion. The ionization 
depth found bv lFortnev et al.l (|2003h is around 1.7 mbar, 
close to the pressure we find for a Na abundance change 
in the condensation model, indicating the effects could be 
important at those pressures and such a model is viable. 

A model with Na ionization effectively lifts the con- 
straint of the temperature profile dropping below the Na 
condensation curve, which the condensation model uses 
to deplete atomic Na. Effectively, a model with Na ion- 
ization leaves a wide range of temperatures possible in 
the middle atmosphere (corresponding in our model to 
T m ) as the collisional broadened Na profile at those al- 
titudes (1,000 - 2,000 km) contains no temperature in- 
formation. The corresponding H2 Rayleigh scattering 
slope (below ~3,400 A) is at the edge of the spectrum 
and also provides no useful middle atmospheric temper- 
ature constraints. Within the framework of equilibrium 
chemistry, the temperature T m is likely below that of 
the TiO/VO condensation curve and above that of the 
Na2S condensation curve. However, Rayleigh scattering 
still implies a high temperature of 2,200±260 K at 33±5 
mbar and the deep narrow Na cores suggest high altitude 
temperatures of around 1,770^570 K. Without a strong 
constraint on the middle atmosphere temperature, T-P 
profile fits are dominated by these upper and lower at- 
mospheric constrains, leading to anomalously high tem- 
perature T-P profiles fit to be everywhere above 2,000 
K. Though technically good fits to the data, such pro- 
files can be ruled out as the temperatures would be ev- 
erywhere higher th an the dayside tempe r ature inversion 
seen w ith Spitzer l|Knutson et al.1 [2008; Bu rrows et al.1 
I2007al) . Additionally, the lack of observed Ca I at A4,227 
A in the low resolution data further suggests tempera- 
tures reach values lower tha n that required for Ca-T i and 
Ca-Al-bearing condensates ()Loders fc Fe glcv 2006j). 

Though the current data lacks enough constraints for a 
unique T-P profile solution in the photoionization model, 
the ionization pressure, Pi on , can be estimated given the 
likely range of temperatures. We fit various isothermal 
models fitting Pi on as free parameter along with the Na 
abundances (taken to be constant) above and below that 
pressure. The resulting ionization pressure varies from 
0.5 mbar to 5 mbar for temperatures of 1,100 and 2,000 
K respectively. Higher temperature fits result in smaller 
upper-atmospheric Na abundances and larger Na abun- 
dance contrasts. The abundances found in the lower tem- 
perature fits are in general agreement with those found 
in the Na condensation model, while the upper temper- 
ature fits find much smaller upper-abundances, with the 
2,000 K model giving an upper Na abundance 180 times 
less than solar. Thus, a Na abundance drop by a factor 
of ten can be considered as a lower limit, and may be 
much larger depending on the depletion mechanism. 

3.3. High altitude temperature inversion 



A ^0.11% absorption of Na indicates that our trans- 
mission spectra probes an altitude range of ~3,500 km. 
Such a large altitude range with our slant transit ge- 
ometry indicates we are sensitive to pressures between 
10's of mbar and ~0.001 mbar. These l imb pressures 
are in agreement with those found from IFortnev et al. 
(2003), who analyzed the original [Charbonneau et al 
( 20021) medium resolution Na measurements. At pres- 
sures higher than 30-50 mbar, our transmission spectra 
is optically thick to Rayleigh scattering by the abundant 
molecule H2, providing a firm upper pressure limit. Such 
pressure restrictions, in addition to the strong narrow Na 
line cores, facilitate the necessity of a hot high-altitude 
temperature (above ~1,500), which then naturally pro- 
vides the observed altitude range due to an increased 
scale height. Our high altitude high limb temperature 
may point toward the atmospheric escape mechanism 
being felt, which necessitates large temperat ure rises to 
perhaps 10,000 or 15,000 K in t he exosphere (|Yelldl2004t 
iVidal-Madiar et al.l[2fM 12001 . 

The highest altitudes probed here are dependent on 
the Na core, however, and thus have a correspond- 
ingly larger associated error. The large error coupled 
with the possible Na depletion mechanisms leaves open 
a wide range of possible upper atmospheric tempera- 
ture profiles. However, a high altitude temperature rise 
can be seen when analyzing the narrowest Na wave- 
length bands alone, sensitive to high altitudes well above 
the condensation/ionization altitude, making the high- 
altitude temperature measurement independent of the 
depletion mechanism. The hot high-altitude tempera- 
ture is illustrated in Fig. 4, which shows pho t omet - 
ric Na measurement s from ICharbonneau et all (2002) 
and ISing et all (|2008[ ) overplotted with various constant 
abundance isothermal models. A single isothermal model 
fits through the 12 and 17 A bands, indicating that the 
temperature and Na abundance do not change signif- 
icantly between those altitudes (^2,000 to 2,200 km). 
However, the absorptions measured in the 4.4 and 8.8 A 
bands, which probe higher in altitude (^2,700 to 3,700 
km), rise much faster than the isothermal model fitting 
the lower altitudes of the 12 and 17 A bands, indicating 
that warmer temperatures are needed. 

In the ionization Na depletion scenario, there are less 
temperature constraints for the middle atmosphere com- 
pared to the Na condensation model. However, limb 
temperatures everywhere above ~ 1,500 K are unlikely, 
and in contradiction with Spitzer measurements of in- 
frared brightness temperature. T-P profiles applicable 
to Na ionization models have less severe temperature 
rises compared to the condensation model, though a hot 
high-altitude temperature is still necessary and at least 
a modest temperature rise probable. 

A similar high altitude temperature rise will also be 
likely necessary to explain recent measurements of an- 
other hot- Jupiter, HD189733b. The Na core in that 
hot- Jupiter reac hes ~4,950 km above the average optica l 
transit altitude (|Redfield et alJ l2008t IPont et al.ll2007h . 
As the atmospheric scale height of HD189733b is smaller 
than HD209458b at the same temperature, a more ex- 
treme high altitude temperature rise is likely necessary, 
however, to fully explain the existing measurements. 
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3.4. Comparison with Spitzer anti-transit data 

A stratospheric temperature inversion has been used 
to interpret Spitzer infr ared secondary eclips e measure- 
ments of HD209458b (jKnutson et alj [2008), thought 
to be caused by an absorber a t optical wavelengths 
on the day-s i de of the planet (|Burrows et al.l [2007a; 
iFortnev et al.l I2007T ) . This optical absorber has 
been postulated to be clouds or molecules TiO and 
VO (iHubenv, Burrows fe Sudarskvl 12003; Bu rrows et al.1 
l2007at IFortnev et al. |[2007|), with a low albedo ruling out 
some species, though not TiO or VO (|Rowe et al.1 120061 : 
iBurrows et"atll2007a[) . 

Compa ring our limb T-P profi le to the day-side pro- 
file from IBurrows et all (|2007aD (see Fig. 3), the at- 
mospheric pressures probed are similar, between 0.01 
and 30 mbar, des pite our slant transit viewing geome- 
try (|Fortne v 2005). The atmosphere in the infrared has 
a higher opacity, due t o strong molecular tra nsitions from 
species such as water ([Burrows et al.ll2007a|) , placing the 
infrared surface of the planet viewed at normal angles at 
similar altitudes as those probed in slant geometry in the 
optical at the limb. 

The high altitude temperature inversion we observe at 
the limb below pressures of ~1 mbar is separate from 
that of the dayside and occurs at a much higher altitude. 
However, the altitude of TiO and VO we find from both 
models, largely confined to pressures ab ove ~10 mbar, 
closely matches the model presented by IBurrows et al.1 
(j2007ah . who placed the stratospheric absorber at pres- 
sures of ^25 mbar and imposed heat redistribution be- 
tween pressures of 10 and 100 mbar. The probable pres- 
ence of TiO and VO in our spectrum along with the alti- 
tude similarity implies TiO and VO are indeed likely to 
be the stratospheric day-sid e temperature inversion insti- 
gators (|Fortnev et al.ll2007f) . As our fit of TiO and VO is 
based largely on low resolution data, a more robust de- 
tection and interpretation of TiO / VO will likely require 
higher resolution HST spectra than those used here, such 
that multiple specific molecular bandheads could be ob- 
served. 

Some form of heat redistribution at or near our low- 
altitude pressure (~30 mbar) would likely be needed to 
explain the high temperature we find there, as the tem- 
peratures there are signif icantly above that pr edicted 
by radiative equilibrium (Showman "eTaLl l20Qgft . Ini- 
tial Spitzer mid-infrared measurements of the bright- 
ness temperature as a function of orbital phase have re- 
vealed that at least 30% of the energy from the stel- 
lar flux must be redistributed to the night side to ex- 
plain th eir null results, implying efficient h eat redistri- 
bution (|Cowan. Agol. fc Charbonneaull2007l ). Compar- 
ing with models by IBurrows et al. I (|2007al lbD with higher 
P n s (>0.35), night side temperatures in the region of 
imposed heat distribution (corresponding to our pres- 
sure at z s ) can have larger temperatures than the day 
side, perhaps explaining our high temperature at z s . Al- 
ternatively, systematic pressure scale differences between 
the two studies may indicate that our transmission spec- 
tra instead probes the top of the stratospheric tempera- 
ture inversion, giving a potential explanation of the high 
temperature at z s as those altitudes could already be 
warm. The minimum temperature we find, however, is 
in general agreement with radiative equilibrium, perhaps 



implying the approximation is valid at those pressures 
around 1 mbar. Such scenarios can be further distin- 
guished in future studies, with tighter constraints on 
the average infrared and optical opacities, additional in- 
frared observations, along with improved global circula- 
tion models. The presence of a temperature inversion, 
TiO and VO along with heat redistribution would seem 
to currently place HD209458b in a intermediate hot- 
Ju piter class betw een the pM and pL classes as proposed 
bv IFortnev et all (|2007l ). 

4. DISCUSSION 

The atmospheric temperatures we find, within the con- 
densation scenario, and presence of TiO and VO can also 
explain other HD209458b observational features. For in- 
stance, the low middle atmospheric temperature and the 
presence of TiO and VO would also explain the lack of 
observed potassium (K), originally predicted to be a ma- 
jor contributor to the optical opacity. At the lower tem- 
peratures of the middle atmosphere, K condenses into 
potassium chloride allowing KC1 to become a dominant 
K-bearing gas. Condensed K would then naturally de- 
plete the atmosphere of atomic K, as seen here with Na, 
leading to reduced signatures. In the lower warmer at- 
mosphere, where wide atomic K line wings would be ob- 
servable even at low resolution, TiO and VO likely make 
up the surrounding continuum further masking the sig- 
nature. If Na is instead ionized, the minimum limb tem- 
perature is higher and ionization of K, which is easier 
than Na to ionize, likely depletes neutral K. If present 
in the atmosphere, a positive detection of atomic K in 
HD209458b will probably require high resolution ground- 
based or medium resolution space-based data to detect 
the K line core. Further detections or upper-limits of 
other atomic species which have different ionization po- 
tentials and condensation temperatures can offer clues 
in future studies when distinguishing between the two 
depletion scenarios. 

The Na mixing ratio we find for the lower atmosphere 
is about two times solar while the abundance in the mid- 
dle and upper atmosphere is found to be 0.2 times solar 
or less, thus our models have decreased Na abundance 
without the need for a global deficiency of Na. In chem- 
ical equilibrium, Na2S condensation effectively removes 
all the Na below the condensation temperature, since 
the protosolar Na abundance is ~13% that of S, limit- 
ing the amount of Na2S which can condense directly to 
the amount of Na. Our condensation model T-P pro- 
file crosses the condensation curve again at higher alti- 
tude, making the altitudes of our minimum temperature 
a "cold trap" for Na. It is currently unclear how the 2 nd 
crossing of the condensation curve could affect the Na 
abundance, as the S/N limits us to fit only a single value 
for the whole upper Na abundance. 

Assuming equilibrium chemistry, our temperature pro- 
file would also result in less CO in the middle atmosphere, 
as our temperature profile at that level is closer to the CO 
= CH4 curve which would result in an enhancement of 
CH 4 and H 2 over that of CO and H 2 (see Fig. 2). This 
is consistent wi th the lack of CO de tection in HD209458b 
during transit ([Deming et al.ll2005h . although deviations 
from equilibrium have been shown to push the CH4/CO 
ratio to sm all values for the pressures p robed here of less 
than 1 bar lCooper fc Showman! (|2006f ). 
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Atmospheric data resulting from transit measurements 
represent the average atmospheric properties across the 
observed limb of the planet, which would include the 
morning, eveni ng, and pole terminato r s. As previously 
recognized by llro. Bezard. fc Guillotl (|2005l ). tempera- 
ture and structure differences could result in differing 
transit signatures between the morning and evening ter- 
minators. Longitudinal differences might also play a 
role, though the transit geometry generally limits sig- 
nificant opacity contributions to regions of the termi- 
nator itself, as the density drops rapidl y from the ter - 
minator along the transit line of sight (Hansen 2008). 
A T-P profile decreasing in temperature and pressure 
with increasing altitude would further mitigate longitudi- 
nal differences, corresponding to pressures above f mbar 
in our condensation model fit. Efficient he at redistri- 
bution (|Cowan, Agol. fc Charbonneau! 120071 ) implies all 
these differences could be minimal, though future Spitzer 
observations and the next generation of global climate 
models would be able to better quantify the differences. 
Given the size of our temperature uncertainties, limb dif- 
ferences of even several hundred degrees would have a 
negligible impact on our resulting model fit. 

Primary transit data alone is currently not able to give 
detailed knowledge on atmospheric differences between 
the morning and evening terminators. Such transit sig- 
natures would appear during ingress and egress, where 
the transit light curve is less sensitive to changes in plane- 
tary radii and more sensitive to errors in limb-darkening. 
Nevertheless, as our Na D line cores remain relatively 
narrow in the medium resolution data, sensitive to the 
middle and upper atmospheres, a solar abundance of Na 
on the potentially warmer evening terminator in the mid- 
dle and upper atmospheres is ruled out. For the conden- 
sation model, this would indicate that Na settles from the 
middle atmosphere and remains at lower altitudes, pre- 
sumably at the Na2S condensation altitude. Condensed 
Na2S would have little impact on anti-transit observa- 
tions and the obs ervational constraint of a low albedo 
(|Rowe et al.ll2006l ). On the day side, the hotter temper- 
atures would sublimate Na2S back into their atomic con- 
stituents, and are likely transpa rent in normal viewing 
geometry anyway (|Fortnevll2005f ). 

5. CONCLUSIONS 

As seen here, fitting primary transit transmission spec- 
tra can give detailed atmospheric information over large 
altitude ranges, and thus large pressure regimes, where 
different physical processes may be important. In our 
case, strong sodium absorption coupled with a relatively 
clear atmosphere provides this large altitude range and 
reveals the pressure where Na depletion via condensa- 
tion or ionization takes place. The depletion level natu- 
rally provides a contrast level between the warmer deeper 
atmosphere, likely sensitive to global circulation, and a 
cooler middle atmosphere where many atomic and molec- 
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ular species have condensed out and/or are depleted via 
ionization. As strong opacity in the core of the Na line 
permits measurements up to very low pressures, allow- 
ing Na to be sensitive to altitudes perhaps as high as 
the lower thermosphere. The slant transit geometry and 
high associated column density restricts the lowest alti- 
tudes probed, however, limiting this study to altitudes 
at and above the stratosphere. 

Interpreting transit spectra becomes easier with in- 
creased wavelength coverage and increased resolution. 
For HD209458b, low resolution data spanning the optical 
reveals a planetary continuum comprised of H2 Rayleigh 
scattering, sodium absorption, and TiO/VO giving a 
self consistent picture of a hot atmosphere at pressures 
around ~10 mbar. This result is consistent with infrared 
Spitzer emission features and a temperature inversion at 
those altitudes. The medium resolution data places the 
detection of Na absorption on firm ground, as there is 
sufficient S /N and resolution to both resolve the Na dou- 
blet and observe the depletion level. This overall atmo- 
spheric picture is not possible without data from all three 
spectroscopic gratings and highlights the importance of 
future observations at other wavelengths and resolutions. 

Transit transmission spectra have become a powerful 
method to probe the abundances of atomic and molecular 
sp ecies for hot-Jupite rs, as seen here for HD209458b and 
in lSwain et all (|2008ft for HDI89733b. Both Jupiter and 
Saturn are enhanced in metals relative to the sun, as 
would seem to be the case here for HD209458b having a 
sodium abundance about twice that of it's parent star. 
Such studies can begin to answer the question of whether 
or not all hot-Jupiters are enriched in heavy elements, 
shedding light on giant planet formation. 

After the next servicing mission, HST will be capable of 
extending the transmission spectra of HD209458b from 
the UV to NIR, potentially linking important physical 
atmospheric processes such as the altitude dependence 
of UV/NUV energy deposition as well as abundance and 
altitudes of other atomic and molecular species. As cur- 
rently a patchwork of such observations exist for the two 
ideal targets HD209458b and HDI89733b, future obser- 
vations should focus on covering wavelength and resolu- 
tion gaps, which can lead to both theoretical and obser- 
vational links between the separate studies. 
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Fig. 1. — Transit Na D absorption from;Sing ct al. (2008). Plotted is the planetary absorption percentage compared to the minimum 
absorption level referenced at 5,000A (1.440%). In black is the low resolution data binned over 16 pixels and overplotted in red is the 
medium resolution data binned over 18 pixels (10 A). Representativ e error bars for the b inned data are also plotted (also see Fig. 10 
k, 11 of Sing et al. 20Ql) as well as the mean planetary radius from Knutson ct al. (2007) (dotted line). Plotted in blue is our transit 
absorption model which contains opacity from atomic Na, Rayleigh scattering, and contains upper-atmospheric Na depletion. The model in 
purple includes Na, TiO, and VO illustrating the wavelengths where those features are prominent. Higher absorption values indicate higher 
at mospheric altitudes, w ith the right y-axis labeled as the altitude above our minimum absorption depth, assuming the system parameters 
of IKnutson et al.l (t200W) . The Na D line cores can easily be seen in the zoomed plot of Fig. 2. The narrow Na Dl and D2 cores extend 
above a wide plateau of Na absorption indicating a low amount of Na is present in the middle atmosphere compared to the lower. The 
drastic abundance change, along with the pressure and temperature where the change occurs, indicates that Na is condensing into sodium 
sulfide on the night side of the planet. 
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Fig. 2. — Same as in Fig. 1, but zoomed on the Na region with the medium resolution data binned over 3 pixels (1.7 A). The Na D 
doublet is easily resolved and reaches up to an average altitude of ~3,500 km. 
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Fig. 3. — Limb temperature-pressure atmospheric profile of HD209458b for the sodium condensation mo del. Plotted (th i ck line , three 
points with error bars) is our best fit terminator T-P profile, while (thin) is a day-side inverted profile from Burrows ct al. (2007a). Also 
shown are the chemical equilibrium condensation curves of iron, TiO, and sodium (dashed lines). The horizontal line represents the pressure 
as which the T-P profile crosses the sodium condensation curve. Above the condensation curve, at high altitudes and low pressures, Na is 
depleted compared to lower altitudes. The CH4/CO line represents the profile where equal amounts of CO and CH4 are expected assuming 
chemical equilibrium. The whole data set with low and medium resolution observations are fitted with a total of 7 free parameters, resulting 
in the plotted T-P profile and sodium abundance below and above the sodium condensation temperature. A fit imposing an adiabatic lower 
T-P profile is also shown. The sodium abundance above the Na2S condensation level is found to be significantly lower than the sodium 
abundance below the condensation level (the ratio of both abundances is significantly above 1 at 3-cr level); this reveals the presence of 
Na2S condensation. 
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Fig. 4. — Differential spectre-photometric measurements of atmospheric sodium absorption from (diamonds) Charbonncau ct al. (2002) 
and (circles) Sing et al. (2008) using the HST/STIS medium resolution data. Also plotted are the Na absorption predictions from isothermal 
models with temperatures (bottom to top) of 500, 750, 1000, 1250, 1500, and 2000 K and constant sodium abundance. The isothermal 
models illustrate the necessity for a high-altitude temperature increase in the narrowest two photometric bands (4.4 and 8.9 A) which probe 
the highest altitudes. 



